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ESTIMATING PROPAGATION VELOCITY THROUGH 
A SURFACE ACOUSTIC WAVE SENSOR 

CROSS-REFERENCE TO RELATED APPLICATIONS 

5 This application claims priority to U.S. Provisional Patent Application 

No. 60/533,177, filed December 30, 2003. 

TECHNICAL FIELD 

The invention relates to surface acoustic wave (SAW) sensors and, more 
10 particularly, techniques for analyzing and interpreting the output of a SAW sensor. 

BACKGROUND 

Chemical and biological testing is commonly used to test for the presence or 
absence of chemical or biological agents. Testing for the presence of chemical or 

15 biological agents in blood, food or other materials is often performed to ensure safety 

or to facilitate diagnosis of medical conditions. For example, testing is used to identify 
chemicals, bacteria or other agents in blood samples taken from medical patients, 
laboratory samples developed for experimental purposes, food samples, or the like. In 
addition, chemical and biological testing is used to test for medical conditions such as 

20 pregnancy, diabetes, bacterial infection, and a wide variety of other conditions that may 

affect the patient's chemistry or biology. 

One type of sensor that has been developed for chemical or biological sensing 
capabilities is a surface acoustic wave (SAW) sensor. One example of a SAW sensor is 
a Love mode shear-horizontal surface acoustic wave (SH-SAW) sensor. A SH-SAW 

25 sensor includes four main components: 1) a piezoelectric stibstrate; 2) an input inter- 

digitated transducer (IDT) on the substrate, which is used to excite an acoustic wave 
based on the piezoelectric effect; 3) an output IDT on the substrate, which receives the 
transmitted acoustic wave and generates electrical output by exploiting the piezoelectric 
effect; and 4) a wave-guide layer over the IDT's, which converts SH-type waves into 

30 waveguide Love modes for transmission from the input IDT to the output IDT. The 

presence of one or more materials on the surface of the SH-SAW sensor affects wave 
propagation through the waveguide layer in response to the presence of bacteria or 
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other agents on the surface of the sensor, which facilitates detection of bacteria or other 
agents. 

SUMMARY 

5 In general, techniques are described for estimating the propagation velocity, or 

equivalently, for estimating a time delay through a surface acoustic wave sensor. In 
particular, techniques which measure and exploit a proper segpient of phase frequency 
response of the surface acoustic wave sensor are described for use as a basis of 
bacterial detection by the sensor. 

10 In one embodiment, the invention provides a method comprising identifying a 

segment of phase frequency response of a surface acoustic wave sensor, and estimating 
a time delay associated with wave propagation through the surface acoustic wave 
sensor based on the identified frequency response. 

In another embodiment, the invention provides a computer-readable comprising 

15 instructions that when executed in a processor identify a segment of phase frequency 

response of a surface acoustic wave sensor, and estimate a time delay associated with 
wave propagation through the surface acoustic wave sensor based on the identified 
frequency response. 

In another embodiment, the invention provides a system comprising a surface 

20 acoustic wave sensor, a sensor analyzer to receive output of the surface acoustic wave 

sensor, and a processor to receive input from the sensor analyzer, identify a segment of 
phase frequency response of a surface acoustic wave sensor, and estimate a time delay 
associated with wave propagation through the surface acoustic wave sensor based on 
the identified segment of phase frequency response. 

25 The invention may be capable of providing one or more advantages. In 

particular, use of change of propagation velocity of a surface acoustic wave sensor, as 
described herein, can improve detection of bacteria via the sensor, relative to 
conventional techniques that use phase shift as the basis for detection. Moreover, use 
of an estimated propagation velocity of a surface acoustic wave sensor as the basis for 

30 detection may allow for detection of bacterial concentration. 

The details of one or more embodiments of the invention are set forth in the 
accompanying drawings and the description below. Other features, objects, and 
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advantages of the invention will be apparent from the description and drawings, and 
from the claims. 

BRIEF DESCRIPTION OF DRAWINGS 

5 FIG. 1 is a perspective view illustrating an exemplary surface acoustic wave 

(SAW) sensor that may be used in one or more embodiments of the invention. 

FIG. 2 is a block diagram illustrating a system according to an embodiment of 
the invention. 

FIG 3 is a flow diagram illustrating a technique according to an embodiment of 
10 the invention. 

FIG. 4 is a graph illustrating an example phase frequency response of a SAW 

sensor. 

FIG. 5 is a depiction of an exemplary LabView screen associated with 
execution of the techniques described herein. 
1 5 FIG. 6-3 1 are various graphs illustrating techniques that can be used in 

accordance with the invention, and various desirable characteristics that can be 
observed relative to the characteristics observed by conventional techniques. 

DETAILED DESCRIPTION 

20 FIG. 1 is a perspective view illustrating an exemplary SAW sensor 10 that may 

be used in accordance with an embodiment of the invention. SAW sensor 10 may 
comprise any of a wide variety of SAW sensors. SH-SAW sensors are typically 
constructed from a piezoelectric material with a crystal-cut and orientation that allows 
the wave propagation to be rotated to a shear horizontal mode, i.e., parallel to the plane 

25 defined by the waveguide, resulting in reduced acoustic damping loss to a liquid in 

contact with the detection surface. Shear horizontal acoustic waves may include, e.g., 
thickness shear modes (TSM), acoustic plate modes (APM), surface skimming bulk 
waves (SSBW), Love-waves, leaky acoustic waves (LSAW), and Bleustein-Gulyaev 
(BG) waves. 

30 In particular, Love wave sensors may include a substrate supporting a SH wave 

mode such as SSBW of ST quartz or the leaky wave of 36°YXLiTaC>3. These modes 
may preferably be converted into a Love-wave mode by application of thin acoustic 
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guiding layer or waveguide. These waves are frequency dependent and can be 
generated if the shear wave velocity of the waveguide layer is dower than that of the 
piezoelectric substrate. 

In one example, sensor 10 comprises a Love mode shear-horizontal surface 
5 acoustic wave (SH-SAW) sensor. 

SAW sensor 10 includes a substrate 12 which typically comprises a 
piezoelectric material. SAW sensor 10 also includes an input inter-digital transducer 
(IDT) 14 on substrate 12, which is used to excite an acoustic wave based on the 
piezoelectric effect. In addition, SAW sensor 10 includes an output IDT 16 on 

10 substrate 12, which receives the transmitted acoustic wave and generates electrical 

output by exploiting the piezoelectric effect. A wave-guide layer 18 is formed over the 
IDT's 14, 16. Wave-guide layer 18 converts SH-type waves into waveguide Love 
modes for transmission from input IDT 14 to output IDT 16. 

A layer of material, such as a layer of antibodies is applied over waveguide 

15 layer 1 8. In operation, a fluid being tested for the presence of bacteria is brought into 

contact with waveguide layer 18. If bacteria is present in the fluid, the bacteria attaches 
to the antibodies on waveguide layer 18 and thereby affect wave propagation through 
waveguide layer 18. Accordingly, analysis of the wave propagation through waveguide 
layer 18 of SAW sensor 10 allows for bacterial detection or detection of other agents 

20 that may interact with a material coated on waveguide layer 1 8. 

In some cases, SAW sensor 10 includes a plurality of sets of input and output 
IDT's. For example, SAW sensor 10 may include a first input IDT 14 and a second 
input IDT 15 which respectively correspond to first output IDT 16 and second output 
IDT 17. In that case, first input and output IDT's 14, 16 comprise the active portion of 

25 SAW sensor 10, and second input and output IDT's 15,17 comprise the reference 

portion of SAW sensor 10. Different types of antibodies may be coated on the surface 
of waveguide layer 18 between first input and output IDT's 14, 16 and second input and 
output IDT's 15, 17, such that only the bacteria of interest bonds to the antibodies 
between the active portion corresponding to first input and output IDT's 14, 16. In that 

30 case, the reference portion corresponding to second input and output IDT's 15, 17 

allows for reference measurements which can account for temperature variance affects, 
or the like, which could otherwise affect wave propagation through SAW sensor 10. 
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SAW sensors are commonly used for bacterial detection but may be designed 
for detection of any of a wide variety of other chemical or biological agents. 
Accordingly, different materials may be coated on the waveguide layer of SAW sensor 
10 in order to facilitate detection of various chemical or biological agents. Waveguide 
5 materials may preferably be materials that exhibit one or more of the following 

properties: low acoustic losses, low electrical conductivity, robustness and stability in 
water and aqueous solutions, relatively low acoustic velocities; hydrophobicity, higher 
molecular weights, highly cross-linked, etc. In one example, SiC>2 has been used as an 
acoustic waveguide layer on a quartz substrate. Examples of other thermoplastic and 
10 crosslinked polymeric waveguide materials include, e.g., epoxy, 

polymethylmethacrylate, phenolic resin (e.g., NOVALAC), polyimide, polystyrene, 
etc. 

In general, the presence of a particular material on the surface of the SAW 
sensor affects wave propagation through the waveguide layer in response to the 

1 5 presence or absence of another material passing over the waveguide layer, which 

facilitates detection of the material passing over the waveguide layer. Accordingly, 
materials coated on the waveguide layer may be selected to attract, trap, bond with or 
otherwise attach to materials suspended in a fluid that flows across waveguide layer 18. 
In this manner, SAW sensor 10 facilitates detection of bacteria or other biological or 

20 chemical agents. 

In order to generate Love mode surface wave, the shear wave velocity of 
waveguide layer 18 is typically lower than that of substrate 12. In that case, the 
acoustic energy will be trapped to near the sensing surface of SAW sensor 10. Love 
mode SH-SAW sensors, in particular, generally have high sensitivity to surface 

25 perturbations. Mass loading can change the surface conditions of a Love mode SH- 

SAW sensor. Therefore, measuring the change of the propagation velocity or resonant 
frequency can be used to quantitatively detect mass loading, and thereby detect the 
presence of a chemical or biological agent. 

One of the most important performance indices for an SH-SAW sensor used as 

30 a detector is the mass sensitivity, i.e., how sensitive is the sensor to the loading mass on 

its surface. For analyzing the mass sensitivity, two indices have been established: 
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the propagation velocities and resonant frequencies of the sensor without and with 
5 surface perturbation resulting from an infinitesimal thin rigid loading layer with mass 

Am s =p m e, where p m and s are the density and thickness of the loading layer. 

SJn » — S f m , where v g is the group velocity, and for Love mode SH-SAW 

> 1 , where the notation vo(o)) as well as v g (a)) is used to emphasize that the 
v g (o)) 8V 

velocities are frequency dependent in the dispersive cases. This means that for Love 

Av 

1 0 mode SH-SAW sensor, using — for detecting the anomalies of the boundary 

v ° 

conditions of the surface of the sensor should be better than using ^ . Accordingly, in 

f o 

Av 

accordance with the invention, — may be used as a basis of a detection indicator for 

SAW-sensor 10, or the like. 

Conventionally, for a Love mode SH-SAW sensor with a triple transit echo 
1 5 (TTE), the phase frequency response was given by: 

or equivalently: 

where f and <o represent frequency and angular frequency, L is the distance between the 
20 center of the input and center of the output IDT, p=a 2 and a is the reflection coefficient 

of the input and output IDTs. Since P is small and can be determined by frequency 
response, without loss of generality, p is assumed to be zero. Thus, for phase frequency 
response, one only needs to be concerned with the case that does not include TTE. The 
phase frequency response is therefore: 
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However, e^ f > = 0-*"™ ^-WW, wherek=0,d;l,±2,+3 etc. Thus, the 
phase frequency response <p(f) is multi-valued and cannot be uniquely determined by 
the value of e J(p(,) . This problem is called phase ambiguity. Only a so-called "main 
5 value" of cp(f) can be determined from e Mf) . Since 

exp{-j2iift(f )} = exp{-j2n(ft(f ) - [fx(f )])} , 
where [x] is the integer part of x, i.e. [x] is the maximum of the integers, which are less 
or equal to x, the main value of the phase response becomes: 

<Mf ) = -2n(fr(f ) - [fi(f )]) + ji (in radian) or 4>(f ) = -360(fi(f ) - [f*(f )j) + 180 (in 

10 degree). 

Thus, -n«|>(f)£n (in radian) or - 1 80<<|>(f)<l 80 (in degree). 
For <p(f) = -2nf-, 



20 



^lv=v 0 = 2*f-^. 

v 0 

Therefore, an approximation of the phase change is given by: 



15 Am = 9*r*_t_ 



Aq> = 2nf— Av. 
Therefore, 
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m = _ 2jl( f t _ [f k }) + „ . if [f J^j = [f — L_j , then £t = -*L . Otherwise, 
v v v 0 v 0 + Av «(> 0 v 0 



-L Av + [f ] - (f — 1 

Aj vg Vp+Av v 0 ^ 

♦o __L + [f _L, *"v 0 



v 0 v 0 



WO 2005/066622 



-8- 



PCT/US2004/042793 



The condition [f— ] = [f — ^_ j implies f— - 



v 0 v 0 +Av 



v 0 v 0 + Av 



<1. Let f|-t--- 



v 0 v 0 + Av 



= 1, 



i.e. f 



LAv 



v 0 (v 0 +Av) 



= 1 . For Av < 0 , it follows that: 



^ _JLandAv = -^- 



v 0 (v 0 +Av) fL fL + v 0 

This means that when lAvl £ — ^ — , [1—] * [f — . For example, If 

fL + v 0 v 0 v 0 +Av 

5 f=l 03MH,v 0 =4000 m/s and L=8.8mm, then [f -^-] * [f Vq ^ Av 1 for |Av) > 1 7.5747 m/s. 

Therefore, when [f — ] * [f — - — ] , A* or — may not be a proper indicator function. 
v 0 v 0 + Av 4>o 

For this reason, the invention provides techniques for estimating the 

propagation velocity through a SAW sensor based on an identified segment of phase 

frequency response. The propagation velocity is simply related to the time delay 

10 through the SAW sensor, as described below. In accordance with the invention, a 

method may include identifying a proper segment of phase frequency response of a 

surface acoustic wave sensor, and estimating a time delay associated with wave 

propagation through the surface acoustic wave sensor based on the identified frequency 

response. 

15 FIG. 2 is a block diagram illustrating a system 20 that includes a SAW sensor 

10, and a sensor analyzer 21 that obtains measurements from SAW sensor 10. System 
20 also includes a processor 22 that interprets the output of SAW sensor 10. In other 
words, sensor analyzer 21 receives output from SAW sensor 10 and provides input to 
processor 22 so that the output of SAW sensor 10 can be interpreted. 

20 Processor 22 receives input from sensor analyzer 21, which comprises 

measurements associated with wave propagation through SAW sensor 10. Processor 
22 then determines whether SAW sensor 10 has detected the presence of particular 
bacteria or other material for which SAW sensor 10 is designed to detect. Processor 22 
executes instructions to perform various techniques and functions ascribed to the 

25 processor herein. Although the invention is not limited in this respect, SAW sensor 10 

may be housed in a cartridge, or the like, and may be electrically coupled to sensor 
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analyzer 21 via insertion of the cartridge into a slot. Processor 22 may be housed in the 
same unit as sensor analyzer 21 or may be part of a separate unit or separate computer. 

Processor 22 may also be coupled to memory 24, which stores a propagation 
velocity routine 25 consistent with the teaching of this disclosure. Alternatively, 
5 propagation velocity routine 25 may be implemented by hardware within processor 22. 

In any case, processor 22 executes propagation velocity routine 25 in order to estimate 
a time delay associated with wave propagating through the surface acoustic wave 
sensor based on the identified frequency response, as described herein. 

By way of example, processor 22 may comprise a general -purpose 

10 microprocessor that executes software stored in memory 24. In that case, processor 22 

may be internally housed in a specifically designed computer, a general purpose 
personal computer, workstation, handheld computer, laptop computer, or the like. 
Alternatively, processor 22 may comprise an application specific integrated circuit 
(ASIC) or other specifically designed processor. In any case, processor 22 executes 

1 5 propagation velocity routine 25 in order to estimate a time delay associated with wave 

propagation through the surface acoustic wave sensor 10 based on the identified 
frequency response, as described herein. 

Memory 24 is one example of a computer readable medium that stores 
processor executable software instructions applied by processor 22. By way of 

20 example, memory 24 may comprise random access memory (RAM), read-only memory 

(ROM), non-volatile random access memory (NVRAM), electrically erasable 
programmable read-only memory (EEPROM), flash memory, or the like. Propagation 
velocity routine 25 such as one of those mathematically described below, are stored in 
memory 24 and may form part of a larger software program used for analysis of the 

25 output of SAW sensor 10. For example, propagation velocity routine 25 may be a sub- 

routine programmed within a Lab View software platform, which is described in greater 
detail below. 

FIG 3 is a flow diagram illustrating a technique according to an embodiment of 
the invention. As shown in FIG .3, fluid is brought into contact with a surface of SAW 
30 sensor 10 (31). The fluid may include samples of material for which bacterial testing is 

needed. For example, the fluid may be allowed to pass over the surface of waveguide 
layer 18 of SAW sensor 10 (FIG 1), which includes antibodies that react to the bacteria 
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of interest. SAW sensor 10, for example, may be housed in a cartridge that defines a 
fluid path over the surface of waveguide layer 18 such that a fluid can be introduced 
into the cartridge and allowed to pass over waveguide layer 18- via the fluid path. 

Processor 22 receives measurements taken from SAW sensor 10 by sensor 
5 analyzer 21, and applies propagation velocity routine 25 stored in memory 24. In doing 

so, processor 22 identifies a segment of phase frequency response of SAW sensor 10 
(32), and estimates a time delay associated with wave propagation through SAW sensor 
1 0 based on the identified frequency response (33). In some embodiments, an 
estimated propagation velocity derived from the estimated time delay can be used to 
10 identify the concentration of the bacteria or other material in the fluid (34). 

Numerous mathematical techniques for estimating the time delay associated 
with wave propagation through the surface acoustic wave sensor will now be discussed. 

T = — refers to the time delay between the center of the input and the center of the 
v 

output IDT separated by distance L. In the dispersive case, the time delay is also a 
15 function of the frequency. That is, x(f) = . It follows from the definition of <|>(f) (in 



degree) that: 



fx(f) = [fT(f)]- 1^ + 0.5. 



For any given running frequency f q (also referred to as the operating 
frequency), there are two phase inflection frequencies f 1 and f 2 proximate to a running 
20 frequency f 0 such that: 

(1) <f 0 <f 2 

(2) <{>(fi) = 180 , *(f 2 -0) = -180 and <t>(f 2 ) = Wi + 0) - 1*° > 

(3) Mfi) = ftxft)], f 2 t(f 2 ) = [f 2 T(f 2 )] and [f 2 x(f 2 )] = [fit(fi)]+l, and 

(4) For any f t <£ f < f 2 ,[fT(f)] = [frff,)] , 

25 FIG. 4 is a graph illustrating an example frequency response of a SAW sensor. 

Labeled in FIG. 4 is an exemplary a running frequency f 0 , a typical first phase 
inflection frequency f-j and a typical second phase inflection frequency f 2 . Again, for 
any fi ^ f < f 2 , [ft(f )] - [fjtft )] . Therefore: 
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Mf) = [fix(f 1 )]-^ + 0.5, forf 1 Sf<f 2 . 
Taking the derivative with respect to f for both sides of the above equality, results in: 

That is, 

which is one basic differential equation for calculating the time delay x(f). 

If the initial condition is: x(f) | f=foo = T(f 00 ) , where f 00 e [f 1f f 2 ) and arbitrary, then 

the solution of the basic differential equation for calculating the time delay x(f) 
becomes: 

10 T (f, = Y^oo)-3^^(f) + 3^^(foo)- 

The solution can be easily derived from: 



and 



^(0 = 1^)1-^ + 0.5 



foot(foo) = [W)]-^r+o.5. 

dOt) 



15 The basic differential equation may be used to prove the following properties of x(f). 

Proposition 1: If for f e [fn.f^) c [f 1f f 2 ) , <K9 is linear with respect to f, then 
x(f) is constant and 

Proof of proposition 1 : 
20 Suppose f 00 e [f 11t f22) • = —AzWoo) is the solution of the basic differential 

360 

equation: 

f *l£> + x(f ) = --L*$yi 9 i n [f 11rf22 ) because the left side 

df 360 df 



WO 2005/066622 



-12- 



PCT/US2004/042793 



,« + T( , ) = t( , ) = -_J-^ 00) and the right side -_J_^l..-J-^ -0 ) due 
to the linearity of <|>(f) in [fn»f22) . It is also from the linearity of <j>(0 inlfn.f^) that 

f 22 _T 11 

Proposition 2: Suppose <()(f) is differentiable in (fi ,f 2 ) and right differentiate 
5 for f-| . Then: 

T(f) = — !_f!*<0 is the first order estimation of x(f). 
360 df 

Proof or proposition 2: 

In the neighborhood of f 00 <= [f 1t f 2 ) , 

♦(f)«Woo) + *(foo)(^-foo). 
10 Therefore, in this neighborhood <j>(f) is approximately linear with respect to f and 

♦(f) - Woo ) • From the proposition 1, x(f ) « -^^*r • T* 1 * is > *< f > = is ^ 

first order of estimation of x(f). 

Algorithms or routines for estimating the Time Delay x(f) will now be 
discussed. From proposition 2, the following algorithms can be used as part of 
15 propagation velocity routine 25 to obtain the estimate of x(f). 

fOOOlJ Algorithm 1: t„ ^and xft) = -3^1) - ~ "X-"" "tt 
Using the equation: 



one obtains: 



1 7 f ff 2 -f, 360 f f 

In particular: 

-„ v_il_J 1 Wo) , 0 5 

( 0, "'o *2-'i 360 f 0 + f 0 ' 

where f<> is running frequency and is also denoted by fnj n in this description. 



WO 2005/066622 



PCT/US2004/042793 

-13- 



10 



Algorithm 2: f 00 = f* , e.g. $(f.) = ±90 . Around f» take 10 frequencies 
f»,f».....ff (10) ^ a linear regression for (f <1) ,«f (1) )).(f< 2 >.ij>(f (2) )),...,(f (10) .«f (10 >)). 

1 

The slope of this line is denoted by<j>(f.)andt(f*) = i(f*) . Therefore, 

360 

W 360 f Vl 1 360 f W # 360 f Yl 1 

Algorithm 3: f 00 = f 0 = 'run • Around fb take 10 frequencies f (1 \f (2) ,... f f (10) 
and a linear regression for (f m .«f (1) )Mf W ,W w )) f . .,(f (10) . Mf m )). The slope of this 
line is 

$(f 0 ) and the first order estimation of x(f 0 ) is i(f 0 ) = -^rW f o) • 

360 

Algorithm 4: Since 



holds for any f 00 e [f 1f f 2 ) - Regarding f 00 as a variable and taking integral with respect 
to foo for both sides of the above approximate equation, provides 

fl 360 f f<| 360 fl f 360 fl f 

which is equivalent to 

15 x(f)(f 2 -fi)--^y(f 2 *f2>-WM - f i )+ ^]f7* (f ^ )d, oo 



i.e., 



Because Wj) = 180 , 4>(f 2 -0) = -180 , 

1 f 1 1 1x/« 0.3 111 f ?. /# ... 
In particular, 
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T(f0> = Uf^h - 360 * + V + 180 U r 00,df00 - 

Comparing the T(f)in algorithm 1 with that of algorithm 4, shows that algorithm 

4 adds one more term ^-^ 1 - / <Mf 00 )df 0 o • If the phase function <J>(f) is symmetric 

180 f f 2 -fi fi 

about the point ((fi+f 2 )/2, 0), then f /<Mf 0 o)dfoo = 0 . In that case, algorithms 1 and 4 are 

5 the same. This means that the term -j-j , 1 - Woo ) df oo is used to reduce the 

180 f f 2 -fi fl 

estimation error of the algorithm 1 resulting from non-symmetry of the phase response 

For the real measured phase response, <|>(f) is only given at the sampling points 
of the frequency f. That is, <K0 is a discrete sequence rather than a continuous function 

10 in the interval [f 1 1 f 2 ) . Furthermore, ) and <t>(f 2 ) may be not equal to 1 80 and -1 80, 

respectively, although it is true that 4>(f-| - Af ) < 0 , <Mf t ) > 0 , $(f 2 ) < 0 and ${f 2 + Af ) > 0 , 
where Af is the sampling interval in the frequency domain. In order to increase the 
estimation accuracy, Af should be as small as possible and linear regression can be used 
to extrapolate <|>(f) from [f l9 f 2 )to [f new1l f new2 ) such that 4>(f ne wi) = 180 311(1 

15 <J>(f n ew2 - 0) = -180 . For the extrapolated <|>(f), 

v 1 f new1 1 1 ±t* x 0-5 1 1 Af f ™2 

*0 f new2 ~ f new1 360 f 0 ?o 180 *0 T now2 -T new1 ( Mw1 

or simply 

n f o) = 7— f } — 7-*ao) + -p- 

f 0 T new2 - T new1 360 T 0 T 0 

Notably, the above algorithms for estimating the propagation velocity are only 
20 based on a proper segment of the measured frequency response of the SH-S AW sensor. 
In at least this way, the algorithms are different than conventional methods. 
Conventional methods, may measure the time delay of the sensor or time delays 
between several fingers of IDT by some time domain device. Alternatively, 
conventional methods may be based on an entire phase frequency response, an 
25 amplitude frequency response, and an inverse Fourier transform. In contrast, the 

techniques described herein are based only on a segment of the measured frequency 
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response of the SH-SAW sensor, i.e., the proper segment, which can be identified as 
described herein. • 

Algorithm 1 was implemented in the software platform Lab View for on-line 
velocity measurement. Lab View is a commercially available software program that can 
5 be obtained from National Instruments, Inc., USA. Using a SAW sensor as a detector 

is generally based on the comparison of propagation characteristics of the sensor 
without and with a particular surface perturbation. Such a comparison can be carried 
out dynamically or statically. 

For the dynamic comparison, the propagation characteristic of the sensor is 

10 measured as a time series. For the static comparison, the sensor should consist of two 

channels, the reference channel and detection channel (also called the active channel). 
For example, the reference channel may correspond to the reference portion between 
IDTs 15 and 17 (FIG. 1), whereas the detection channel may correspond to the active 
portion between IDTs 14 and 16. The surface perturbation takes place in the detection 

1 5 channel. In that case, the propagation characteristics of two channels of the sensor are 

measured. The measurement of propagation characteristic of SH-SAW sensor is the 
base for using the sensor as a detector. 

Resonance frequency frcsonancc and propagation velocity v are two of the most 
direct propagation characteristics, and Av/v are more sensitive than Af/f to the surface 

20 perturbation. Therefore, the measurement of the propagation velocity may provide an 
improved method for efficiently using SH-SAW sensor as a detector. In contrast, 
conventional techniques typically measure the log-amplitude response A(f) and phase 
response <j>(f) of SH-SAW sensor by a sensor analyzer and read out the phase at the 
running frequency from the measured phase frequency response. Developing routines 

25 for the direct measurement of the propagation velocity of the SH-SAW sensor becomes 

very important in the practical application of SH-SAW sensor. The algorithms 
presented above provide a software-based approach to estimate the propagation 
velocity from a proper segment of the measured phase frequency response, e.g., which 
can be measured by an 8753ET network analyzer commercially available from Agilent 

30 Technologies, Inc., USA. Other similar network analyzers could also be used. 

LabView was chosen as any easy way to both communicate with the sensor 
analyzer and to do the calculation based on the model. An exemplary LabView screen 
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associated with execution of the techniques described herein is shown in FIG. 5. The 
inputs are in part defined by the spectral resolution of the sensor analyzer and include 
the minimum and maximum frequency (equal to a total of 1600 data points on the 
sensor analyzer) and the running frequency. 
5 The outputs from Lab View allow various interpretations of the results to be 

made. For the majority of the plots shown on the Lab View screen shown there are two 
outputs, one representing the measurement sensor and the other representing the 
reference sensor. In Lab View, plots of phase and velocity can be used to determine if a 
sensor is working reliably. If the phase becomes noisy or if the amplitude decreases 
10 more than about 20db then the sensor is deemed bad and discarded. 

The phase histogram plot can be used as another measure of sensor quality and 
relates to non-linearity in the phase plot. The velocity plot may provide the result of 
the estimation according to the techniques described herein. 

Plots of magnitude and phase are taken at the running frequency and plotted 
15 against time. Lastly, temperature and analog values are outputs that are plotted against 

time. In order to calculate velocity from a proper segment of phase frequency response 
an internal sub- VI is used within Lab View. The first step is to calculate all the +180 
and -180 degree phase points. In other words, to determine the phase inflection 
frequencies, the technique may sample a plurality of phase responses at frequencies 
20 proximate to the running frequency and estimate the phase inflection frequencies as a 

function of the plurality of phase responses. 

The difference between the velocity expression in continuous frequency and its 
discrete realization was discussed above. However, this first step may also be applied 
during the implementation of Algorithm 1 in Lab View. For example, if the sensor 
25 analyzer has a limited resolution in frequency, the data about a -180 to 180 degree 

phase transition may be similar to that below: 

-175.97, -176.866, -177.761, -178.657, -179.552, 179-5522, 178.6567, 
177.7612, 176.8657, 175.9701, 

An extrapolated linear fit may be used to find the exact -1 80 and 1 80 degree 
30 phase points. An illustration of exemplary + 1 80 and -1 80 degree phase points is 

shown in FIG. 6. FIG. 7 illustrates use of a bilaterally extrapolated linear fit to find the 
exact -180 and 180 degree phase. As can be appreciated from FIG. 7 it is possible to 
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calculate the values of fi, fi and the phase at the running frequency <(>(fb). However, 
because the format is different in Lab View it may be necessary to use different 
nomenclature from the model derived above. For example, in Lab View: 

Tau_Fr = (((Fl/(F2-Fl))*360)-Pha - Fr+180)/F_run*1060000*360) (with 
5 F_run converted to hertz), 

and then, 

V_Fr=(L*0.001)/Tau_Fr (with L converted from millimeters to meters). 
The final output of velocity (V_Fr) is plotted in Lab View and corresponds to v(f 0 ) , 
discussed above. 

Using a separate application of Lab View, it is possible to take stored phase files 
and reconstruct velocity data. The phase data can then be transformed using exactly the 
same model as above. In this way, it is possible to see the effect of changing the 
running frequency relative to the velocity obtained. Lab View allows two running 
frequencies (at 0.1MHz difference) to be evaluated simultaneously. Advantageously, 
the phase inflection frequencies define edges of a monotonically changing subset of a 
graph of phase versus frequency of the surface acoustic wave sensor. Thus, the 
techniques described herein may allow not only for the identification of a material via 
the sensor, but also an indication of the concentration of the material. 

Example 1 

The Love mode SH-SAW sensor used in the experiment was a LiTaCh device 
operating at 103 MH Z provided by Sandia National Laboratories, USA. A low- walled 
flow cell was placed over the sensor and filled with Phosphate Buffered Saline (PBS) 
buffer solution at pH 7.5. This liquid container was connected to a syringe pump 
system to allow a slow flow of buffer. During the experiment, multiple aliquots of 250 
microliters of 0.05 mg/ml Bovine Serum Albumin %(BSA) protein was injected into 
the cell at designated times. An 8753ET network analyzer from Agilent Technologies, 
Inc., USA, measured the log-amplitude and phase response of the sensor approximately 
every twenty seconds. Based on the resulting phase responses, the curve of phase at the 
operation frequency versus time was immediately obtained. Algorithm 1 was then used 
to calculate the propagation velocity versus time at the operating frequency. FIGS. 8 
and 9 respectively show the phase response and propagation velocity plotted over the 
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time of the experiment. The valleys of the curve of velocity versus time precisely 
correspond to moments when BSA was injected into the very slowly flowing buffer 
stream. Compared with the curve of phase versus time, the calculation of propagation 
velocity provides much more information about the mass loading on the surface of a 
5 Love mode SH-SAW sensor, resulting in an increase in sensitivity to mass loading and 

viscosity. 

Example 2 

In this example, three separate concentrations of bacteria being detected were 
10 used. The first bacterial concentration was 10 3 cfu/ml (colony forming unit per 

milliliter), and the results are illustrated in FIGS. 10 and 11. The second bacterial 
concentration was 10 5 cfu/ml and the results are illustrated in FIGS. 12 and 13. The 
third bacterial concentration was 10 7 cfu/ml and the results are illustrated in FIGS. 14 
and 15. 

Ay 

15 In particular, the graphs of FIGS. 10, 12 and 14 show— versus time at the 

specified concentrations, and the graphs of FIGs. 11,13 and 15 show ^ versus time 

s>o 

respectively at the specified concentrations. The usual way to use an SH-SAW sensor 
for mass-loading measurement is that the sensor contains two channels: reference 
channel and detection channel (the detection channel also being referred to herein as 

20 the "active channel" or "active portion" of the sensor). In the reference channel, the 

surface of the sensor is largely unperturbed, while it is perturbed in the detection 
channel. By measuring and comparing the propagation characteristics of the two 
channels, the perturbation resulting from the mass-loading in the detection channel can 
be detected. If the indicator function, which can be used to describe the difference of 

25 the propagation characteristic of the two channels, is monotonic with respect to the 

mass of the loading material, then the existence of loading mass can be detected and 

Ay 

quantified. It has been pointed out above that using — as indicator function for 

detection is better than using — . This example demonstrates the correctness of this 
conclusion from a real experiment. In this experiment, the surface of the detection 
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channel was provided with bacterial loading at specific concentration. The 
concentrations of the bacteria used in the experiment are 10 3 , 10 5 , and 10 7 cfu/mL, 
respectively. The delay from injecting bacteria into the detection channel until back- 
washing was around 20 minutes. The curves of — and versus time for each 

5 bacterial concentration are shown in FIGS. 10-15, where the calculations of the 

propagation velocities of the two channels are based on the above algorithm. It should 
be pointed out that the experiment in this example includes no* antibody. Bacteria were 
bound directly to the surface chemistry which was a saccharin immobilization layer. 

The — and ^ at the steady state are extracted from the above curves for each 
v 0 4>o 

10 bacterial concentration. The relationships of — versus log-concentration and 

v o 

^ versus log-concentration are respectively shown in FIGS. 16 and 17. 

Although the experimental results are just for three different bacterial 
concentrations, clearly, — is not monotonic with respect to the bacterial concentration. 

It also demonstrates the correctness of the conclusion that ^ is not a proper indicator 

<Po 

1 5 for quantitatively detecting loading mass with large dynamic range. 



Example 3 

The Love mode SH-S AW sensor used in the experiment was a LiTa0 3 device 
operating at 103 MH Z provided by Sandia National Laboratories, USA. During the 
20 experiment, the surface of the sensor was: 

1 . In air for approximately 20 min 

2. In buffer for approximately 30 min 

3. Subjected to bacteria in the buffer for approximately 18 min 

4. Briefly washed and left in buffer for approximately 20 min 
25 5. Subjected to De Ionized (DI) water for approximately 30 min. 

Bacteria were bound to the antibody attached to the surface. FIG. 18 is a graph 
of the propagation velocity at the operating frequency versus time and FIG. 19 is a 
graph of phase at the operating frequency versus time. The propagation velocities are 
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calculated using algorithm 1 presented above. In FIGS. 18 and 19, the dashed vertical 
lines define the time durations of various surface conditions of the sensor as detailed 
above. 

By observing the graph of velocity versus time in FIG. 18, the process and 
5 results of the experiment can be understood much more easily. The propagation of the 

surface acoustic wave of the sensor has higher velocity when its surface was in air as 
shown at 181. Upon injecting buffer, the velocity decreased and tended to a steady 
state as shown at 182. Upon adding a certain amount of bacteria, the velocity 
decreased again as shown at 183. The brief washing did not remove the bacteria, which 

10 were adsorbed by the antibody attached to the surface of the sensor. Therefore, the 
velocity during washing as that in the previous state as shown at 184. Finally, the 
velocity of the surface acoustic wave sensor increased because the surface condition 
was changed from the adsorbed bacteria to water as shown at 185. By stark contrast, it 
is difficult to understand and explain the experimental process and results based on the 

15 curve of phase versus time shown in FIG. 19. 



Example 4 

The Love mode SH-SAW sensor used in the experiment was a LiTa03 device 
operating at 103 MH Z provided by Sandia National Laboratories, USA. A low-walled 

20 flow cell was placed over the sensor and filled with Phosphate Buffered Saline (PBS) 

buffer solution at pH of approximately 7.5. This liquid cgntainer was connected to a 
syringe pump system to allow a slow flow of buffer. During the experiment, multiple 
aliquots of 250 microliters of Bovine Serum Albumin (BSA) protein at various 
concentrations were injected into the cell at designated times. The injection of proteins 

25 conformed to Table 1, below. In Table 1, BSA refers to Bovine Serum Albumin, and 

IgG refers to Immunoglobulin. 
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TABLE 1 



Injection 

moment 

(sampling 

points in 

time 

domain) 


Injected 
Proteins 


Concentration 
(mg/ml) 


Age of Solution • 


Locations of 
jumping 
down point 
of velocity 


Av/v 


/u 


no A 

doA 


a ac 


1 day old 


/ I 


A Al O"* 


OA 


IgVJ 




O w raal/O aI/1 

z weeKs oia 






1 

1ZZ 


TIC A 


a 


l aay oia 


1 OA 
1ZO 


A A1 OA 


1 <o 


r>uner 










1 QA 


"DC A 


a a< 


2 weeks old 






Tin 


TJG A 


A A< 


1 aay old 


HA 1 
/41 


A A1 OO 


/oJ 


TJC A 


A AO< 


1 day old 


to / 


A m AO 
-U.UIOZ 


Q 1 A 
olO 


TJC A 

JboA 


A A1 


i aay oia 


OZU 


n nun 

-U.U1 1U 




HQ A 


A A17** 
U.Ujjj 


l aay oia 


OJO 


n ni ao 

-U.U IU7 


ooo 


TJC A 


A A/L 


i aay oia 


QQO 


A A1 OA 
-U.U 1 oo 


916 

/1U 




0 05 


1 dav aIH 


920 


-0.0201 


948 


BSA 


0.00625 


1 day old 


952 


-0.0056 


994 


BSA 


0.003125 


1 day old 


997 


-0.0011 


1023 


BSA 


0.05 


1 day old 


1027 


-0.0201 


1058 


IgG 


0.05 


1 day old 


1061 


-0.0014 


1084 


BSA 


0.05 


1 day old 


1087 


-0.0209 



As described in Example 1, an 8753ET network analyzer from Agilent 
Technologies, Inc., USA, measured the log-amplitude and phase frequency response of 
5 the sensor approximately every 20 seconds. Based on proper segments of these phase 

frequency responses for the active channel, and algorithm 1 proposed above, both the 
curves of phase and propagation velocity at the running frequency versus time fare 
were obtained. FIG. 20 is the resultant graph of phase versus time, whereas FIG. 21 is 
the resultant graph of propagation velocity versus time. 
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The details of the curve of propagation velocity (m/s) versus time (sampling 
point in time domain) for BSA injections are shown in the Figures 22-24. The locations 
of the valleys of the curve are very close to the injection moments listed in Table 1. In 
FIG. 23 the labels correspond to the following: 
5 (1) (737) BSA .05 mg/ml; 

(2) (783) BSA .025 mg/ml; 

(3) (816) BSA .0125 mg/ml; 

(4) (854) BSA .0333 mg/ml; 

(5) (888) BSA .04 mg/ml; 
10 (6) (916) BSA .05 mg/ml; 

(7) (948) BSA 0.00625 mg/ml; 

(8) (994) BSA 0.003125 mg/ml; 

(9) (1023) BSA .05 mg/ml; 

In order to estimate Av/v more accurately, some segments of the curve of 

15 velocity around the injection moments of BSA are used. An example is shown in FIG. 

26. In this case, the estimated injection moment is the 952 nd sampling point in time 
domain. From this segment illustrated in FIG. 26, the edges of the valley are 
determined and the velocities v 0 and Vt are estimated. Then, the Av/v=(vi-v 0 )/v 0 for 
this injection is calculated, in this case, Av/v=-0.0056. FIG. 25 shows the 

20 corresponding segment of the curve of <J>. It is difficult to determine A<J>/<|> for this 

injection from the graph of FIG. 25. 

More examples are presented in FIGS. 27-30. In particular, FIGS. 28 and 30 
illustrate segments of the curve of v around the estimated injection moments and FIGS. 
27 and 29 illustrate the corresponding segments of the curve of <|>. The calculated 

25 (Av/v)'s for FIGS. 28 and 30 are -0.0201 and -0.0197, respectively. It is still difficult 

to determine A<J>/<j> from FIGS. 27 and 29. 

The more accurate estimations of Av/v about the estimated injection moments 
of BSA are presented in Table 1. The resulting Av/v versus the known concentration of 
the injected BSA at the same injection moment, which is listed in Table 1, is shown in 

30 FIG. 31 as a circle. Corresponding to the concentration 0.05 mg/ml, there are seven 
overlapping (Av/v)'s. 
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. One can examine the circle corresponding to the largest Av/v. The curve of 
FIG. 3 1 results from connecting the circles. The curve of FIG. 3 1 is monotonic. 
Moreover, the monotonic nature of the curve of FIG. 31 demonstrates that using Av/v, 
which can be determined by the curve of velocity and calibration, the sensor may 
further detect the concentration of the injected BSA if the concentration of the injected 
BSA is in the interval [0.003125, 0.05] mg/ml. 

The measurement of propagation characteristic of SH-SAW sensor is the base 
for using the sensor as a detector. The resonance frequency fiance and propagation 
velocity v are two direct propagation characteristics and Av/v is more sensitive than 
A£7f to the surface perturbation. Therefore, the measurement of the propagation 
velocity is a preferred index for efficiently using SH-SAW sensor as a detector. 

In contrast, the conventional approach is to measure the log-amplitude 
frequency response A(f) and phase frequency response <|>(f) of SH-SAW sensor by a 
sensor analyzer, and then read out the phase at the running frequency. The serious 
limitation of A<j> or A<|>/<|> for detection was discussed above. By developing hardware 
and/or software for the direct measurement of the propagation velocity of the SH-SAW 
sensor, a practical application of a SH-SAW sensor can be achieved. Based on the 
system model of SH-SAW sensor with conventional triple transit echo (TTE) 
techniques, several algorithms have been presented to estimate the propagation velocity 
from a proper segment of the phase frequency response measured by a sensor analyzer. 
The first algorithm (algorithm 1) has also been implemented in the software platform 
Lab View for on-line velocity measurement. 

Compared with the conventional phase-based or phase shift- based method, the 
techniques described herein will not lose any information carried by phase. Also, 
because of overcoming the phase ambiguity, the techniques described herein can 
extract much more information about the surface condition of SH-SAW sensor relative 
to conventional techniques. Furthermore, the examples showed the monotonic nature 
of the curve of Av/v versus concentration of the liquid mass loading in a large dynamic 
range of the concentration. This may therefore be the basis for using SH-SAW sensor 
in practice to quantitatively detect mass loading. In contrast, for the curve of A<J>/<j> 
versus concentration of the liquid mass loading, monotonic curves only hold true for 
small concentrations. 
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Various embodiments of the invention have been described. In particular, 
techniques for estimating the propagation velocity through a surface acoustic wave 
sensor have been described. The techniques may be implemented in hardware, 
software, firmware, or the like. Example hardware implementations include 
5 implementations within a general purpose microprocessor, an application specific 

integrated circuit (ASIC), a field programmable gate array (FPGA), specifically 
designed hardware components, or any combination thereof. In addition, one or more 
of the techniques described herein may be partially or wholly executed in software. In 
that case, a computer-readable medium may store or otherwise comprise computer- 

10 readable instructions, i.e., program code, that can be executed.by a processor or to carry 

out one of more of the techniques described above. The techniques can be used for 
both constant time delay (in non-dispersive case) and frequency dependant time delay 
(in dispersive case). 

For example, the computer-readable medium may comprise random access 

1 5 memory (RAM), read-only memory (ROM), non-volatile random access memory 

(NVRAM), electrically erasable programmable read-only memory (EEPROM), flash 
memory, or the like. These and other embodiments are within the scope of the 
following claims. 

Suitable methods for coating the devices of the present invention include 
20 Applicants' Copending Application U.S. Serial No. 10/607,698, filed June 27, 2003. 

The present invention may be utilized in combination with various materials, 
methods, systems, apparatus, etc. as described in various U.S. patent applications 
identified below, all of which are incorporated by reference in their respective 
entireties. They include: U.S. Patent Application Serial Nos. 60/533,162, filed on 
25 December 30, 2003 ; 60/533, 1 78, filed on December 30, 2003; 1 0/896,392, filed July 

22, 2004; 10/713,174, filed November 14, 2003; 10/987,522, filed November 12, 2004; 
10/714,053, filed November 14, 2003; 10/987,075, filed November 12, 2004; 
60/533,171, filed December 30, 2003; 10/960,491, filed October 7, 2004; 60/533,176, 

filed December 30, 2003; 60/533,169, filed December 30, 2003; , titled 

30 "Method of Enhancing Signal Detection of Cell- Wall Components of Cells", filed on 

even date herewith (Attorney Docket No. 59467US002); , titled "Soluble 

Polymers as Amine Capture Agents and Methods", filed on even date herewith 
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(Attorney Docket No. 59995US002); , titled "Multifunctional Amine 

Capture Agents", filed on even date herewith (Attorney Docket No. 59996US002); 

PCT Application No. , titled "Surface Acoustic Wave Sensor Assemblies", 

filed on even date herewith (Attorney Docket No. 58928WO003); PCT Application No. 

5 , titled "Acousto-Mechanical Detection Systems and Methods of Use", 

filed on even date herewith (Attorney Docket No. 59468 WO003); PCT Application No. 

, titled "Detection Cartridges, Modules, Systems and Methods", filed on 

even date herewith (Attorney Docket No. 60342 WO003); and PCT Application No. 
, titled "Acoustic Sensors and Methods", filed on even date herewith 

10 (Attorney Docket No. 60209 WO003). 

The complete disclosures of the patents, patent applications, patent documents, 
and publications cited herein are incorporated by reference in their entirety as if each 
were individually incorporated. Various modifications and alterations to this invention 
will become apparent to those skilled in the art without departing from the scope and 

15 spirit of this invention. It should be understood that this invention is not intended to be 

unduly limited by the illustrative embodiments set forth herein and that such 
embodiments are presented by way of example only, with the scope of the invention 
intended to be limited only by the claims. 



